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TO WHOM IT MAY CONCERN: 

Be it known that Kullervo Hynynen, a Citizen of Finland of Medfield, MA and 
Gregory T. Clement, a U.S. Citizen of Boston, MA have an invention entitled SHEAR 
MODE THERAPEUTIC ULTRASOUND of which the following description in 
connection with the accompanying figures is a specification. 



SHEAR MODE THERAPEUTIC ULTRASOUND 



CROSS-REFERENCE TO RELATED ACTIONS 

This application claims the benefit of U.S. Provisional Application No. 60/463,589 filed 
5 April 17, 2003 and entitled "Transskull Shear Mode Ultrasound" which is incorporated here by 
reference. 

STATEMENT AS TO FEDERALLY-SPONSORED RESEARCH 

This invention was made at least in part with Government support under Grant No. NIH 
10 CA76550, awarded by the National Institutes of Health. The Government has certain rights in 
this invention. 

FIELD OF THE INVENTION 

The invention relates to ultrasound diagnostic and therapeutic applications using shear 

15 waves. 

BACKGROUND OF THE INVENTION 

Since its inception in the early 1980's, transcranial ultrasound Doppler imaging has 
demonstrated an ability to measure blood flow, hemorrhages, and perfusion in the brain. Recent 

20 research has also examined the possibility of diagnosing certain degenerative disorders such as 
Parkinson's and depression. These approaches use transcranial Doppler sonography (TCD), or 
the related transcranial color coded sonography (TCCS) to record frequency shifts in the sent and 
backscattered signals. Signals are generally applied with a relatively low frequency (~2 MHz) 
probe in order to penetrate the skull bone, and are often used in conjunction with a contrast agent. 

25 There seems to be little advantage to applying frequencies above 2 MHz, as the skull's increased 
attenuation at higher frequencies causes the bone to act as a low pass filter, returning only the 
lower spectral frequencies. Center frequencies at and below 1 MHz have also been examined, 
showing a stronger signal intensity, but with an expected reduced resolution. 
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The primary setbacks for transcranial procedures are strong attenuation and the distortion 
caused by irregularities in the skull's shape, density, and sound speed. These properties 
collectively contribute toward destroying an ultrasound focus and/or decreasing the ability to 
spatially register received diagnostic information. 
5 Coherent noninvasive focusing of ultrasound through the human skull has been suggested 

for a number of therapeutic and diagnostic implications in the brain. For example, ultrasound 
has been considered as a tool for the transskull treatment of brain tumors, targeted drug delivery, 
improved thrombolytic stroke treatment, blood flow imaging, detecting internal bleeding, and 
tomographic brain imaging. Although the human skull has been the barrier to clinical realization 

10 of many of these applications, studies have demonstrated both minimally invasive and 
noninvasive aberration correction methods for transskull focusing. Minimally invasive 
approaches use receiving probes designed for catheter insertion into the brain to measure the 
amplitude and phase distortion caused by the skull and then correct the beam using an array of 
ultrasound transducers. Alternatively, a completely noninvasive approach uses X-ray computed 

15 tomography (CT) images to predict the longitudinal wave distortion caused by the skull. 

Noninvasive focusing with a therapeutic array has been demonstrated with a longitudinal wave 
propagation model, but the amplitude of the focus was observed to drop when the focus was 
directed close to the skull surface. 

The assumption that the transcranial propagation is composed of mainly longitudinal 

20 modes is valid for small incident beam angles, but rapidly breaks down beyond approximately 
25°, as the longitudinal wave approaches its critical angle. This is a plausible explanation for 
reduced amplitude using the longitudinal model - As the focus is directed toward the periphery 
of the brain, an increasing number array elements are oriented at higher incident angles to the 
skull. 

25 Modeling of shear waves has been dismissed as being either of insignificant amplitude or, 

if significant, that the resulting beam would be incoherent and hard to predict. The absence of 
significant information on the skull bone's elastic wave speeds also has inhibited its 
consideration in modeling. Similar issues exist with respect to ultrasound propagation through 
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other bony structures. 

SUMMARY OF THE INVENTION 

It has been discovered that ultrasound beams can be effectively, coherently transmitted 
5 through bone via shear waves. This can be accomplished, for example, by increasing the angle 
between the beam and the normal to the bone surface beyond the critical angle (approximately 
20°) that causes complete reflections of the longitudinal waves. At angles between 
approximately 25° and 60°, shear waves are generated and propagated through the bone and at the 
inner surface converted again into longitudinal waves that propagate into the soft tissues on the 

10 other side of the bone. This finding has several implications: First, the shear waves have 

roughly the same propagation velocity in the bone as the longitudinal waves have in soft tissues. 
Thus, the wave front may not be significantly distorted by the variable thickness of bone as it is 
for longitudinal waves due to the higher speed of sound in bone (2500-3000 ms/s). It is possible 
to focus ultrasound beams though bone without performing complex patient-specific aberration 

15 corrections. For example, diagnostic ultrasound imaging of the brain can be performed as well as 
ultrasound brain treatments and diagnosis with relatively simple devices. A second implication 
is that the entrance angle of the beam can be increased and can reach larger brain volumes with 
higher gain than using longitudinal waves. Third, predictions of ultrasound focusing through a 
skull can be significantly improved by combining shear wave propagation modeling with 

20 longitudinal propagation models. Further, undesired reflections from bone layers may be 
reduced or eliminated, making diagnostic information easier to analyze. 

Embodiments of the invention may provide one or more of the following capabilities and 
can be used for one or more of the following uses. The invention can be used in ultrasound 
diagnosis and/or therapy and can aid methods that rely on transmitting ultrasound through bone. 

25 In particular, the invention can be applied to systems that seek to image blood flow through bone, 
or can be an enhancement/add-on for existing ultrasound imaging and/or therapeutic devices, or 
in a stand-alone device for transbone ultrasound therapy and/or applied to techniques that rely on 
opening of the blood-brain barrier (BBB) using ultrasound or any other diagnostic or therapeutic 
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procedures conducted using ultrasound. Focusing of ultrasound could be greatly simplified, and 
phase distortion correction could be reduced and possibly eliminated. The tissue volume that can 
be currently reached with ultrasound can be increased. In addition, embodiments of the invention 
can detect cavities in bony structures and/or detection/determination of cavity contents with 
5 undesired reflections from bone surfaces being reduced or eliminated. 

In general, in an aspect, the invention provides a method of diagnosing a subject by 
delivering ultrasound signals using shear waves, the method including applying a portion of an 
ultrasound mainbeam to a bone surface at an incident angle relative to the surface of the bone to 
induce shear waves in the bone, energy in the shear waves forming a substantial part of energy of 

10 first ultrasound waves at a desired region in the subject through the bone, detecting at least one of 
reflected and scattered energy of the applied ultrasound mainbeam, and analyzing the detected 
energy for a diagnostic purpose. 

Implementations of the invention may include one or more of the following features. The 
portion of the ultrasound mainbeam is applied to the bone surface between a longitudinal wave 

15 critical angle associated with the subject and a shear wave critical angle associated with the 
subject. The analyzing includes producing an image of at least a portion of the desired region. 
The desired region is a linear region along a line of transmission of the ultrasound mainbeam. 
Applying the portion of the ultrasound mainbeam comprises applying the portion of the 
ultrasound mainbeam to bone. The bone is a skull, and the portion of the ultrasound mainbeam 

20 is directed at the skull at the incident angle in order to reach the desired region within the skull. 
The desired region is one of a sinus cavity and an inner ear cavity, the method further comprising 
providing an indication of whether the desired region is at least partially fluid filled. The desired 
region is one of a tooth and a jaw bone, the method further comprising at least one of: providing 
an indication of whether the desired region has at least one of a cavity and an abscess, and 

25 obtaining an image of at least one of anatomy and blood perfusion of the desired region. 

Implementations of the invention may also include one or more of the following features. 
The at least a first ultrasound mainbeam is applied in multiple pulses of different frequencies. 
The different frequencies are within a range of frequencies from about 0.3MHz and about 5MHz. 
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The different frequencies are within a range of frequencies from about 1MHz and about 3MHz. 
The pulses have durations within a range of about 1 cycle to about 100 cycles. The different 
pulses have corresponding different amplitudes to compensate for different attenuation amounts 
associated with the different frequencies. The at least a first ultrasound mainbeam is applied in a 
5 burst of multiple cycles, at least two of the cycles having at least one of different frequencies, 
different phases, and different amplitudes. 

In general in another aspect, the invention provides a system for diagnosing a subject by 
delivering ultrasound signals to a target region in the subject using shear waves, the system 
including a transceiver configured to transmit ultrasound energy, directing means, coupled to the 

10 transceiver, for causing a portion of a mainbeam of the transmitted ultrasound energy to be 
incident upon a bone surface of the subject such that ultrasound energy in the portion of the 
mainbeam from the source will induce shear waves in the bone with energy in the shear waves 
forming a substantial part of energy of ultrasound waves at the target region, and analyzing 
means, coupled to the transceiver, for analyzing energy from the portion of the mainbeam that is 

15 returned from the target region for diagnostic purposes. 

Implementations of the invention may include one or more of the following features. The 
directing means are configured to cause the portion of the mainbeam to be incident upon the 
bone surface at a first angle between a longitudinal critical angle associated with the bone and a 
shear critical angle associated with the bone. The directing means comprises at least one of (1) a 

20 positioner configured to mechanically direct a normal direction associated with the source toward 
the surface of the bone at the first angle, (2) a phase/delay adjuster, where the transceiver 
comprises a plurality of radiating elements, the phase/delay adjuster being configured to regulate 
at least one of phases and delays of the plurality of radiating elements to electronically steer at 
least the first mainbeam, and (3) an actuation regulator, where the transceiver comprises a 

25 plurality of radiating elements, the actuation regulator being configured to actuate the elements at 
different times to direct the mainbeam as desired. The positioner is configured to at least one of 
(1) couple to the subject and the transceiver in a fixed manner such that the normal is directed 
toward the surface at the first angle, and (2) mechanically adjust the transceiver such that the 
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normal is directed toward the surface at the first angle. The transceiver comprises a plurality of 
elements configured to radiate ultrasound energy, the system comprising a controller configured 
and coupled to cause at least a portion of the transceiver to emit ultrasound energy, to process 
indicia of returned energy due to the emitted energy to determine an orientation of at least a 
portion of the surface relative to the transceiver, and to actuate only elements of the source that 
have their mainbeams at least partially directed at the portion of the surface between the 
longitudinal critical angle and the shear wave critical angle. The controller is configured to 
process the indicia of returned energy to form an image of the at least a portion of the surface. 

Implementations of the invention may also include one or more of the following features. 
The system includes a controller coupled to the transceiver and configured to actuate the 
transceiver to produce the mainbeam for transmitting energy to the target region, and the surface 
of the subject is an outer surface of a skull of the subject. The system includes a controller 
coupled to the transceiver and configured to actuate the transceiver to produce the mainbeam in a 
plurality of pulses with different frequencies and amplitudes. The different frequencies are 
within a range of frequencies from about 0.5MHz and about 5MHz. The different frequencies 
are within a range of frequencies from about 1MHz and about 3MHz. The pulses have durations 
within a range of about 1 cycle to about 100 cycles. The different amplitudes compensate for 
different attenuation amounts associated with the different frequencies. The system includes a 
controller coupled to the transceiver and configured to actuate the transceiver to produce the 
mainbeam in a burst of multiple cycles, at least two of the cycles having at least one of different 
frequencies, different phases, and different amplitudes. The analyzing means is configured to 
provide an indication of whether the target region is at least partially fluid filled. 

In general, in another aspect, the invention provides a system for delivering ultrasound 
signals to a target region in a subject using shear waves, the system including a transducer device 
configured to transmit and receive ultrasound energy, a controller coupled to the transducer 
device and configured to actuate the transducer device to transmit ultrasound energy toward the 
subject, and a positioning device coupled to the transducer device and configured to ensure that a 
portion of a first mainbeam from at least a portion of the transducer device is directed at a portion 



of a surface of the bone at an incident angle between normal incidence and a shear critical angle 
associated with the subject such that ultrasound energy in the first mainbeam will induce shear 
waves in the subject and energy from the first mainbeam will reach the target region, with energy 
in the shear waves forming a substantial part of energy of ultrasound waves at the target region, 
5 where the controller is configured to cause the transducer device to transmit energy in a at least 
one of: a single pulse including multiple cycles, at least two of the cycles having at least one of 
different frequencies, different phases, and different amplitudes, and a plurality of pulses, with 
different pulses having a different frequencies and amplitudes, and where the controller is 
configured to analyze ultrasound energy from the first mainbeam that is returned from the target 
10 region and is received by the transducer device to determine diagnostic information from the 
returned energy. 

Implementations of the invention may include one or more of the following features. The 
transducer device comprises a plurality of elements configured to transmit ultrasound energy, and 
the controller is configured to inhibit actuation of at least one of (1) a portion of the transducer 

15 device configured to produce a second mainbeam that would be incident upon the portion of the 
surface of the subject at a second angle that is less than the longitudinal critical angle, and (2) a 
portion of the transducer device configured to produce a third mainbeam that would be incident 
upon the portion of the surface of the subject at a third angle that is greater than the shear wave 
critical angle. The different frequencies are within a range of frequencies from about 0.1MHz 

20 and about 5MHz. The different frequencies are within a range of frequencies from about 

0.31MHz and about 3MHz. The pulses have durations within a range of about 1 cycle to about 
100 cycles. The positioning device is configured to couple to the subject to mechanically orient 
the source relative to at least one of the subject and the bone as desired. The transducer 
comprises a plurality of elements configured to transmit ultrasound energy, and the positioning 

25 device is configured to affect phases of the elements to electronically steer the first mainbeam. 
The transducer comprises a plurality of elements configured to transmit ultrasound energy, and 
the controller is configured to affect timing of actuations of the elements to electronically steer 
the first mainbeam. The controller is configured to provide an indication of whether the target 



region is at least partially fluid filled based upon the diagnostic information determined by the 
controller from the returned energy. 

Implementations of the invention may also include one or more of the following features. 
The controller is configured to cause a second mainbeam to be directed at the target region to 
5 stimulate motion in the target region, and w the controller is configured to provide an indication 
of the motion in the target region based upon the diagnostic information determined by the 
controller from the returned energy. The first and second mainbeams have different frequencies. 
The controller is configured to cause second and third mainbeams to be directed at the target 
region to stimulate motion in the target region, and the controller is configured to provide an 
10 indication of the motion in the target region based upon the diagnostic information determined by 
the controller from the returned energy. The controller is configured to produce an image of at 
least a portion of the target region from the returned energy. The target region is a linear region 
of the subject and the controller is configured to produce a linear image from the returned energy. 
In general, in another aspect, the invention provides a method of delivering ultrasound 
1 5 signals using shear waves, the method including applying a portion of at least a first ultrasound 
beam to a subject at at least a first incident angle relative to the surface of the subject to induce 
shear waves in the subject, energy in the shear waves forming a substantial part of energy of first 
ultrasound waves at a desired region in the subject at a therapeutic level. 

Implementations of the invention may include one or more of the following features. The 
20 portion of the first ultrasound mainbeam is applied to a surface of the subject between a 
longitudinal wave critical angle associated with the subject and a shear wave critical angle 
associated with the subject. Applying the portion of at least a first ultrasound beam comprises 
applying ultrasound energy at multiple incident angles between the longitudinal wave critical 
angle associated with the subject and the shear wave critical angle associated with the subject to 
25 focus ultrasound energy in the desired region. The method further includes applying at least a 
portion of a second ultrasound beam to the subject the subject to induce shear waves in the 
subject and to produce second ultrasound waves in the subject at the desired region, producing an 
image of at least a portion of the desired region, and identifying, from the image, whether 
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ultrasound energy from the portion of the second ultrasound beam reaches the desired region in a 
desired manner. Applying the portion of at least a first ultrasound beam comprises applying the 
portion of at least a first ultrasound beam to bone. The bone is a skull, and the portion of the first 
ultrasound beam is directed at the skull at the at least a first incident angle in order to reach the 
5 desired region within the skull. 

Implementations of the invention may include one or more of the following features. The 
portion of at least a first ultrasound beam is applied in multiple bursts of different frequencies. 
The different frequencies are within a range of frequencies from about 0.1MHz and about 5MHz. 
The different frequencies are within a range of frequencies from about 0.2MHz and about 

10 3MHz. The pulses have durations within a range of about 1 cycle to continuous wave. Applying 
the portion of at least a first ultrasound beam comprises applying a portion of a third ultrasound 
beam to the subject to produce shear waves in the subject to produce third ultrasound shear 
waves in the desired region. The portion of the third ultrasound beam is separate from the 
portion of the first ultrasound beam. Applying the portion of at least a first ultrasound beam 

15 comprises applying a portion of a fourth ultrasound beam to the subject at a fourth incident angle 
that is less than the longitudinal wave critical angle associated with the subject. 

In general, in another aspect, the invention provides a system for delivering ultrasound 
signals to a target region in a subject using shear waves, the system including a source configured 
to transmit ultrasound energy, and directing means, coupled to the source, for causing a portion at 

20 least a first mainbeam of the transmitted ultrasound energy to be incident upon a surface of the 
subject to induce shear waves in the subject, energy in the shear waves forming a substantial part 
of energy of first ultrasound waves at the target region in the subject at a therapeutic level. 

Implementations of the invention may include one or more of the following features. The 
directing means is configured to direct the first mainbeam at the surface of the subject at a first 

25 angle between a longitudinal critical angle associated with the subject and a shear critical angle 
associated with the subject. The directing means comprises at least one of (1) a positioner 
configured to mechanically direct a normal direction associated with the source toward the 
surface at the first angle, and (2) a phase/delay adjuster, where the source comprises a plurality of 
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radiating elements, the phase/delay adjuster being configured to regulate at least one of phases 
and delays of the plurality of radiating elements to electronically steer at least the first mainbeam. 
The positioner is configured to at least one of (1) couple to the subject and the source in a fixed 
manner such that the normal is directed toward the surface at the first angle, and (2) mechanically 
5 adjust the source such that the normal is directed toward the surface at the first angle. The source 
comprises a plurality of elements configured to radiate ultrasound energy, the system comprising 
a controller configured and coupled to cause at least a portion of the source to emit ultrasound 
energy, to process indicia of reflected energy due to the emitted energy to determine an 
orientation of at least a portion of the surface relative to the source, and to actuate only elements 

10 of the source that have their mainbeams at least partially directed at the portion of the surface 
between the longitudinal critical angle and the shear wave critical angle. The controller is 
configured to process the indicia of reflected energy to form an image of the at least a portion of 
the surface. The system includes a controller coupled to the source and configured to actuate the 
source to produce the first mainbeam and a second mainbeam at least a portion of which would 

15 be incident upon a surface of the subject at a second angle between the longitudinal critical angle 
associated with the subject and the shear critical angle associated with the subject such that 
ultrasound energy in the second mainbeam from the source will induce shear waves in the subject 
and energy from the second mainbeam will reach target region, where the second angle is 
different from the first angle. 

20 Implementations of the invention may include one or more of the following features. The 

system includes a controller coupled to the source and configured to actuate the source to 
produce the first mainbeam for transmitting energy to the target region. The system includes a 
controller coupled to the source and configured to actuate the source to produce the first 
mainbeam in a plurality of pulses with different frequencies. The different frequencies are within 

25 a range of frequencies from about 0. 1MHz and about 5MHz. The different frequencies are 

within a range of frequencies from about 0.2MHz and about 3MHz. The pulses have durations 
within a range of about 1 cycle to continuous wave. The system includes a controller coupled to 
the source and configured to actuate the source to produce the first mainbeam and a third 
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mainbeam at least a portion of which would be incident upon a surface of the subject at a third 
angle that is less than the longitudinal critical angle associated with the subject. 

In general, in another aspect, the invention provides a system for delivering ultrasound 
signals to a target region in a subject using shear waves, the system including a source configured 
5 to transmit ultrasound energy, a controller coupled to the source and configured to actuate the 
source to transmit ultrasound energy toward the subject, and a positioning device coupled to the 
source and configured to ensure that a portion of a first mainbeam from at least a portion of the 
source is directed at a portion of a surface of the subject at a first angle between a longitudinal 
critical angle associated with the subject and a shear critical angle associated with the subject 
10 such that ultrasound energy in the first mainbeam will induce shear waves in the subject and 
energy from the transmitted ultrasound will reach the target region, where the controller is 
configured to cause the source to transmit energy in a plurality of pulses, with each pulse having 
a different frequency. 

Implementations of the invention may include one or more of the following features. The 
15 source comprises a plurality of elements configured to transmit ultrasound energy, and the 
controller is configured to inhibit actuation of at least one of (1) a portion of the source 
configured to produce a second mainbeam at least a portion of which would be incident upon the 
portion of the surface of the subject at a second angle that is less than the longitudinal critical 
angle, and (2) a portion of the source configured to produce a third mainbeam at least a portion of 
20 which would be incident upon the portion of the surface of the subject at a third angle that is 
greater than the shear wave critical angle. The different frequencies are within a range of 
frequencies from about 0.1MHz and about 5MHz. The different frequencies are within a range 
of frequencies from about 0.2MHz and about 3MHz. The pulses have durations within a range 
of about 1 cycle to continuous wave. The positioning device is configured to couple to the 
25 subject to mechanically orient the source relative to the subject as desired. The source comprises 
a plurality of elements configured to transmit ultrasound energy, and the positioning device is 
configured to affect phases of the elements to electronically steer the first mainbeam. 

Various aspects of the invention may provide one or more of the following capabilities. 
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Ultrasound can be propagated through bone, e.g., the skull, while experiencing reduced distortion 
and increased signal strength, allowing clearer and more accurate brain images, and better 
focusing for therapy applications compared to previous techniques. Unwanted echoes from 
overlying bone layers may be reduced, making analysis of diagnostic information received from 
5 target tissue easier. Bone cavities and/or their contents may be detected. Transbone shear mode 
• propagation can be used for a number of imaging problems including vessel detection, tumor 
detection, tissue morphology, and hemorrhaging in the brain. Transbone ultrasound can be 
provided with reduced distortion and/or higher location accuracy than previous techniques. 

These and other capabilities of the invention, along with the invention itself, will be more 
10 fully understood after a review of the following figures, detailed description, and claims. 

BRIEF DESCRIPTION OF THE FIGURES 

FIG. 1 is an image of a portion of a skull and a simplified line drawing of layers of the 

skull. 

15 FIG. 2 is a schematic diagram of an ultrasound therapy system according to the invention. 

FIGS. 3-9 are schematic diagrams of exemplary portions of the system shown in FIG. 2 in 
use in various applications. 

FIG. 10 is a block flow diagram of a process of using the system shown in FIG. 2 for 
ultrasound diagnostic applications. 
20 FIGS. 1 1-12 are two different mounting procedures used in experiments for transskull 

measurements. 

FIG. 13 is a pair of graphical representations of experimental and simulated pressure 
amplitude and phase as functions of incident angle. 

FIGS. 14A-B are plots of pressure amplitude and phase as a function of incident angle of 
25 a transskull experiment. 

FIG. 15 are images of ultrasound transmitted through a skull at various incident angles. 

FIG. 16 is a plot of measured and simulated field measurements of ultrasound transmitted 
through a skull at a 32° incident angle. 
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FIGS. 17A-B are images of non-normalized and normalized coded excitations. 

FIGS. 18A-B are plots of summed, filtered and unfiltered, respectively, transducer 
response to excitations shown in FIG. 17. 

FIG. 19 shows A-line scans of shear waves and longitudinal waves. 

FIG. 20 shows images of a nylon nut at 0° and 33° incident angle through skull. 

FIG. 21 is a schematic diagram of an experimental setup for testing cavity imaging using 
shear wave propagation of ultrasound. 

FIGS. 22A-B are plots of responses to shear waves and longitudinal waves, respectively, 
incident upon water-filled and air-filled cavities. 

FIG. 23 is a graphical representation of pressure amplitude due to longitudinal and shear 
waves as functions of incident angle upon a plastic layer. 

FIG. 24 is a graph of an exemplary, narrow mainbeam antenna pattern. 

FIG. 25 is a graph of an exemplary, broad mainbeam antenna pattern. 

FIG. 26 is a block flow diagram of a process of using the system shown in FIG. 2 for 
ultrasound therapeutic applications. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Embodiments of the invention provide techniques for transskull and other transbone 
propagation that deliberately induces a shear mode in bone. For transskull propagation, incident 
waves experience a mode conversion from an incident longitudinal wave into a shear wave in the 
bone layers and then back to a longitudinal wave in the brain. The skull's shear speed may 
provide a better acoustic impedance match, less refraction, and less phase alteration than its 
longitudinal counterpart. Using a shear wave, ultrasound may be focused in the brain. 
Longitudinal waves will not be induced in the bone if the ultrasound is incident upon the bone at 
an angle beyond Snell's critical angle. Numerical analysis is provided and demonstrations of the 
phenomena were studied with plastic phantoms and using an ex vivo human skull. Embodiments 
of the invention may be used for various applications, including therapeutic and diagnostic 
applications, which are discussed separately below. Other embodiments are within the scope of 
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the invention. 

It has been observed that under certain conditions it is possible to propagate ultrasound 
through the skull with reduced distortion and higher signal amplitudes by using high incident 
angles. Both numeric and experimental investigation indicate that this is due to the behavior of 
5 shear modes induced in the skull bone. When the ultrasound angle of entry is beyond Snell's 
critical angle for the longitudinal pressure wave, propagation in the bone is purely due to a shear 
wave. This conversion from a longitudinal wave (skin) to a shear wave (skull) and back to a 
longitudinal wave (brain) does not necessarily produce a highly distorted or small-amplitude 
wave. Investigations that studied isotropic phantom materials and later ex vivo human bone 

10 samples support this hypothesis. Indeed, it has been discovered that a focused beam, traveling as 
a shear wave in the skull, may be less distorted than a longitudinal one. In some instances a 
focused beam was observed to be larger in amplitude than a longitudinal wave propagated 
through the same skull area. The success of the longitudinal-shear-longitudinal propagation may 
be primarily due to the similarities between the elastic (shear) wave speed (-1400 m/s) and the 

1 5 sound speeds of water (-1 500 m/s), skin (~1 525 m/s), and the brain (-1550 m/s) and other soft 
tissues. In contrast, the range of longitudinal sound speeds in the skull at relevant frequencies is 
approximately twice these values. Thus, intentionally producing shear-modes in the skull bone 
may be used as a mechanism for producing brain images, or for detecting abnormalities within 
the brain. Similar techniques may be used for therapeutic and/or diagnostic uses in other parts of 

20 the body, or in applications aside from human or other animal bodies. 

Embodiments of the invention use shear-modes intentionally produced in the skull bone 
as a mechanism for producing or enhancing a focus through the skull. Shear modes are 
incorporated into a transcrainal propagation model, demonstrating an ability to predict ultrasound 
phase and amplitude at high incident angles. To illustrate the approach, the field through a single 

25 plastic layer is first modeled. The approach is then applied to sections of human bone, which are 
compared to experimental measurements. The approach could be used for more accurate 
focusing in the brain and could extend the focusing region beyond current methods. 
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Theoretical Foundation 

Based on the skull's small curvature relative to ultrasound wavelengths, the inner and 
outer surfaces of the skull are sectioned into regions that are each approximated as flat, but that 
are not necessarily parallel. This approximation is made to simplify the analysis and is not 
5 required for implementation of the invention. A diagram of the problem is given in FIG. 1, 
showing the ultrasound field divided into areas incident upon these regions of the skull surface. 
Each region is modeled as a single-layer isotropic solid with Lame constants determined as a 
function of the mean density over the propagation region. 

Propagation through an arbitrarily oriented isotropic skull layer is achieved by spectral 

10 decomposition of the incident wave and determination of the ray paths, and attenuation as a 
function of angular wavenumber as the field crosses the skull. Each of these initially 
longitudinal harmonic spectral wave components will be considered below in terms of its 
velocity potential. Without loss of generalization a given component may be viewed in a 
reference frame where the surface normal is oriented along Cartesian y-axis and the z-axis is 

15 defined by the unit vector of the cross product between the propagation axis and the surface 
vector. In this frame, the velocity potential is expressed as 

, I _ aI i(tot-k[xs\\\§' -kyycosQ 1 ) w f(GW-*'jrsin6 7 +k l y ycosQ i ) 

r ~ A L e + A LR e ? (1) 

where A L is the amplitude of the longitudinal wave component incident upon the layer surface, 
A LR is the amplitude of the reflected longitudinal wave, 6 is the angle of incidence and k x and k y 
20 are the wavevector components in the specified reference frame. The superscripts I to HI are 
used to denote to the skin, skull and brain, respectively while the subscripts L and S refer to 
Longitudinal or Shear waves. Accordingly, the transmitted longitudinal potential in the skull is 
given by: 

<f> =A L e (2) 
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and the shear vector potential is 

?'=AZe n ^* Bg -* 3 ~" !i z. ■ (3) 

Using this description, each spectral wave component must be viewed in its own unique 
reference frame. At the boundary, the incident wave is split into a reflected wave, a transmitted 
5 longitudinal wave and a transmitted shear wave. The amplitudes of these waves may be 
determined relative to the incident wave using methods outlined by Kino (Acoustic Waves: 
Devices, Imaging, and Analog Signal Processing, Englewood Cliffs, New Jersey: Prentice-Hall, 
1987). Specifically, the normal component of the particle displacement, 
r =V<|)+Vx\j7 (4) 

10 must be continuous at the boundary as well as the normal stress 
dr dr 

S„=\ 1 + + (X + 2») 1 ± (5) 

ox oy 



and the shear stress 



— + — 
{dy dxj 



(6) 



with the shear and longitudinal sound speed of a given medium related to the Lame constants p. 
15 and X by 

c = £ 

Y m 

Details of the amplitude calculations are provided in Appendix A. Each plane wave solution of 
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describes the behavior of a single angular wavenumber. These amplitudes are calculated for each 
component wavevector space. However, this may be readily performed in closed form, as given 
in Appendix B. 

After propagating into the skull, the longitudinal and shear waves are treated separately, 
5 with the total wave reaching the brain then equal to 

AlU _ ^///^/((oz-^rsine^-^jvcose 777 ) A m ( "(G>r-*f*sinef -*;Vosef ) 

r ll ^ h ls e 5 ( o) 

where A"[ and A™ s are the amplitudes of the longitudinal waves due to the incident longitudinal 
and shear waves in the skull. 

The values of the velocity potentials may be found by equating Eqs (4)-(6) at the skin- 

10 bone interfaces after substituting in Eqs. (l)-(3), and solving for A" and A" . At the skull-brain 
interface, the incident shear and longitudinal velocity potential amplitudes will be equal to the 
product of these transmission amplitudes and the absorption loss experienced within the bone. 
Since each spectral component will have its own independent path length through the skull, its 
total absorption will generally differ between components. A two dimensional representation of 

1 5 the problem is given in FIG. 1 . 

To find the wave amplitudes in the brain, A ! " L and A® , the problem is once again reduced 
to two dimensions by rotating the problem into a reference frame where the surface normal is 
oriented along Cartesian y-axis and the relevant wavevector lies in the x-y plane. Since the soft 
tissue of the brain is fluid-like, the incident longitudinal wave in the skull bone will be further 

20 divided into a reflected shear wave, a reflected longitudinal wave and a transmitted longitudinal 
wave. The shear wave in the skull will be similarly divided, but with differing reflection and 
transmission angles. 

The acoustic pressure at a given point the brain is obtained by separately solving for 
amplitude and phase of each spectral component over a planar area. The spectrum is inverse- 
25 transformed to give the pressure over the measurement plane. The pressure amplitude is equal to 
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the negative of the normal stress which is related to the amplitude of the scalar velocity potentials 
in the brain, A™ and A% by Eq. (4) and Eq.(5). The pressure amplitudes were calculated as 
described in Appendix A. FIG. 23 shows the angular dependence of A l [[ and A% resulting from 
an infinite plane wave in water after traveling through an ideal plastic (e.g., acrylic) layer. 
5 Based on she skull's high attenuation coefficient, additional contributions due to multiple 

reflections within the skull bone are neglected. The total acoustic pressure at any point in the 
brain can be determined by calculating the path length and total transmission amplitude for each 
spectral component. The relative acoustic phase and the overall attenuation (including absorption 
loss) are calculated. The algorithm for determining the field path lengths from the transducer to 
10 the measurement point in the brain is presented in Appendix B. 

Diagnostic Applications 

Referring to FIG. 2, an ultrasound diagnostic system 10 includes an imager 12, a phased 
array 14 of n transducer elements 16, a signal adjuster 18, a controller 20, a frequency generator 

15 22, and a positioner 23. The system 10 is configured to provide ultrasound diagnostic 

capabilities. The system 10 is configured to determine one or more characteristics of an object, 
here a skull 28 of a patient 30, and to apply ultrasound energy (e.g., in the range from about 0.01 
MHz to about 10 MHz , and preferably between about 0.01 MHz and about 3MHz) that is 
focused within the object, here in the patient's brain. While the discussion here focuses on a 

20 skull, other objects may be used, e.g., areas away from the brain, such as sinus cavities, ear canal, 
etc. Further, the invention may be used in dental applications (teeth), application through the 
ribs, spine, or any other bone. The invention could be used to determine if a cavity in a bone 
contains air or fluid and/or the viscosity of the fluid. In addition, the invention may be used to 
image bone marrow or nerves and vessels in the bones. Imaging of the spinal canal is also 

25 possible. The invention may also be used to detect bone cavities, fractures and/or tumors. The 
invention may also be used to determine properties of bone by comparing the shear wave speed 
and the longitudinal wave speed or by applying a force (e.g., an ultrasound-induced radiation 
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force or a mechanical force) on the bone and using ultrasound to detect displacement. 

Signals to the array 14 are provided by a driving arrangement. This arrangement may be 
similar to that reported in Daum et al., "Design and Evaluation of a Feedback Based Phased 
Array System for Ultrasound Surgery." IEEE Trans. Ultrason. Ferroelectr. Freq. Control 
5 45(2):431-4, 1998, but with a driving frequency selected between about 0.1 MHz and about 10 
MHz. The driving can be also performed by other techniques that provide desired ultrasound 
signals. The power and phase to each transducer element 16 may be manually controlled or 
automatically controlled using software and feedback. While an array of elements is preferred, 
the system 10 may be used with a non-coordinated set of elements, or even a single element, 

1 0 depending upon the application. 

The array 14 of transducer elements 16 is configured to be disposed on or near the 
external surface of the patient's head or other surface overlying a bone 28. The array 14 is 
configured in a curved shape, e.g., spherical, although sections of other shapes are possible such 
as planar, including linear. The array 14 is configured to be placed on or near the patient's head 

15 and for focusing ultrasound energy at a distance from the surface of the array 14. The array 14 
can be of a variety of sizes, e.g., 30cm in diameter or larger to fit around a person's head, or 
about 1mm in diameter or even smaller. The elements 16 are preferably transmitters configured 
to pulse, and direct, ultrasound at power levels that depend on the application, and in addition 
receive ultrasound energy. The transducer elements 16 of the array 14 are preferably 

20 piezoelectric transducer elements arranged in the array 14 as shown. The transducer elements 16 
(e.g., piezoelectric ceramic pieces with a center-to-center spacing of A/2) are mounted in silicone 
rubber or other material suitable for damping mechanical coupling between the elements 16. 
Other materials may also be used for the array construction. For example, the array 14 may be 
formed from one or more pieces of piezocomposite material, or any material that converts 

25 electrical energy to acoustic energy. 

The array 14 is coupled to the signal adjuster 18 that is further coupled to the frequency 
generator 22. The frequency generator 22 is configured to provide a common radio frequency 
(RF) signal as the input signal to the signal adjuster 18. Also a separate signal generator could be 
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used for some or all of the elements 16. The radio frequency generator 22 can be of any type that 
will produce the appropriate signals for the signal adjuster 18. Individual frequencies can also 
used. Using an individual signal generator (for example, digital waveform generators), the 
phase, amplitude, and delay of each signal can be set by the signal generators without an 
5 additional circuit. The excitations cause the ultrasound energy to transmit through the patient's 
skull 28 and, depending upon the relative delay, to focus the energy at a selected region within 
the patient's brain if desired. The generator 22 is coupled to the adjuster 18 to split the 
generator's output signal to provide n input signals to the signal adjuster 18. 

Coupled to receive each of the n input signals from the frequency generator 22 are n pairs 

10 of amplifiers 24i-24 n and associated delay circuits 26 r 26 n of the signal adjuster 18. Each pair of 
delay circuit 26 and amplifier 24 represents a channel of the signal adjuster 18. The delay 
circuits 26 are configured to provide n independent output signals to the amplifiers 24 by altering 
or adjusting the delay (and possibly amplitude) of the incoming signals from the generator 22 by 
respective delay factors <Xi-a n . The amplifiers 24j-24 n are configured to amplify the signals from 

15 the delay circuits 26 and to provide the amplified signals to the transducer elements 16 through 
connections, e.g., coaxial cables, individually connecting the amplifiers 24 and the transducer 
elements 16. 

The delay factors a r a n of the delay circuits 26 provide steering of the ultrasound beam 
absent an object in the path of the ultrasound energy. The phase delay for each element 16 

20 associated with steering can be computed using known techniques. The delays oti-a n are 

provided by the controller 20. The controller 20 includes logic that may be provided by software, 
hardware, firmware, hardwiring, or combinations of any of these. For example, the controller 20 
can be a general purpose, or special purpose, digital data processor programmed with software 
instructions in a conventional manner in order to provide and apply the delays oci-a n to the delay 

25 circuits 26, although other configurations may be used. 

The positioner 23 is configured to help ensure that the array 14 can provide ultrasound 
energy to the object 28 at angles to induce shear waves in the object 28. The positioner 23 may 
be shaped to orient the array 14 (or a single element 16) with respect to the outer surface of the 
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object 28 such that a normal direction relative to the array 14, or single element 16, is between a 
longitudinal critical angle (Snell's angle) plus a lower safety factor and below the shear critical 
angle (above which shear waves will not be propagated into the object 28) minus an upper safety 
factor. For a soft tissue-skull interface, Snell's angle is about 20° and the shear critical angle is 
5 about 65° and thus ultrasound may be preferably directed at a skull between about 25° and about 
60°. The safety factors may help account for a width of the mainbeam of the array 14 or element 
16. Referring to FIGS. 24-25, a mainbeam 200 of the array 14 may be relatively narrow, e.g., 
±3° at 3dB down from the normal at 0°. A mainbeam 202 of a single element 16 may be much 
broader, e.g., ±30° at 3dB down from the normal at 0°. For diagnostic applications, the 

10 mainbeam is preferably narrow whether it is from an array (large or small), or a single element. 
The positioner 23 helps to ensure that at least a desired portion of the mainbeam 200, 202 can be 
directed at the surface of the object 28 at an incident angle between the longitudinal and shear 
critical angles. For example, for a single element 16, the positioner 23 may be a block 
configured to receive the element 16 and to be placed on the surface of the object 28 to align the 

15 element 16 as desired. For example, the positioner 23 could be configured to be placed on a 
person's face adjacent the person's nose and to orient the element 16 such that at least a desired 
portion (e.g., with at least a threshold amount of power) of the mainbeam 202 will be directed 
into the skull of the person at an angle between the longitudinal and shear critical angles (e.g., 
toward the person's sinus cavity). Also, shallower angles between the critical angle for the 

20 longitudinal wave and the normal may be used. 

The controller 20 is configured to direct ultrasound from the array 14 as desired. The 
controller 20 is configured to mechanically and/or electronically steer/direct ultrasound energy 
from the array 14. The controller 20 can send instruction signals to the positioner 23 to 
mechanically alter the position of the array 14 and thus the physical pointing of the array 14 (i.e., 

25 the direction of the mainbeam(s) from the array 14 and/or location(s) at which the mainbeam(s) 
is(are) incident upon the object 28). The controller 20 can further control the phase provided by 
the phase shifters 26, or the timing of excitations of the elements 16, to electronically scan the 
ultrasound beam from the array 14. Several techniques could be used for steering by controlling 
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the excitation timing, e.g., the controller 20 could provide delayed excitation signals, or the 
controller 20 could provide a common excitation signal that is time delayed different amounts for 
different elements 16, etc. The controller 20 is configured to ensure that at least a desired 
amount of energy from the ultrasound mainbeam, be it mechanically and/or electrically steered, 
5 is directed into the object 28 between the longitudinal and shear critical angles. The ultrasound 
energy is converted from longitudinal waves into shear waves and transmitted in the object 28 as 
shear waves. Some energy, however, may be directed at the object 28 at angles above or below 
the range between longitudinal and shear critical angles. Further, the shear waves may be 
converted back to longitudinal waves in the object 28 if there is a transition in the object 28 (e.g., 
10 a change in acoustic impedance) that would induce such a conversion (e.g., a transition from 
skull to brain). 

The mechanical direction/orientation of the array 14 may also be manually adjusted. 
Preferably, the angle of the surface of the object 28 is determined by sight, without the use of 
imaging equipment such as a CT (computer tomography) or MR (magnetic resonance) scanner. 

15 Guides, however, may be provided for indicating when a beam normal to the array 14 will 
produce a shear wave in the object 28. The controller 20 is further configured to selectively 
excite the elements 16. Thus, the controller 20 can selectively excite elements that will provide 
an incident angle to the object 28 that will produce a shear wave in the object 28 and not excite 
elements that will not produce a shear wave in the object 28. Further, multiple beams may be 

20 produced simultaneously from the array 14 directed at different locations, and of similar or 
different frequency and/or amplitude. 

The controller 20 is configured to control the delay circuits 26 and the amplifiers 24, and 
to process data received from the object 28, here the skull 28, to provide diagnostic abilities. The 
controller 20 can cause the emitted ultrasound to be sent along a line within the object 28, and to 

25 be scanned (e.g., moved methodically) about a broader region of interest in the object 28. The 
line through the object 28 taken by the ultrasound may change directions (e.g., at interfaces of 
varying acoustic impedance). The controller 20 can process reflected energy received by the 
array 14 into an image based on the direction of the beam and the time between sending and 
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receiving the energy. Also, the controller 20 does not have to process received energy to form an 
image of the object 28. The controller 20 can process the received energy to provide other 
diagnostic information, e.g., a binary indication of whether a void is detected (e.g., a cavity or 
abscess in a tooth or bone, such as a jaw bone), or the presence or absence of material in a void 
5 (e.g., fluid in a sinus cavity or inner ear), etc. The controller 20 can also process the information 
to form an image, e.g., of anatomy of the object 28, of blood perfusion in the object 28, etc. 

The controller 20 is further configured to control the frequency generator 22, the 
amplifiers 24, and the delay circuits 26 to affect the frequency, magnitude, and phasing of 
ultrasound energy supplied by the array 14 to the object 28 for diagnostic applications. It has 

10 been discovered that a non-uniform excitation scheme helps overcome transmission attenuation, 
thus improving energy delivery and reception characteristics. The controller 20 is configured to 
cause the array to be excited with a series of extended bursts, thus producing a series of extended 
ultrasound bursts transmitted to the object 28. The bursts are preferably repeated over a range of 
frequencies and returned signals numerically combined. For example, 10 cycles of energy can be 

15 sent into the skull 28 at frequencies ranging from about 0.5MHz to about 2.0MHz at intervals of 
about 0.01MHz. The amplifiers 24 are preferably controlled such that desired (e.g., uniform) 
signals are sent (or received) at each component in the spectrum. Also, power of the transmitted 
signals is preferably higher at higher, more attenuating frequencies (e.g., above about 1MHz for 
skulls), for example being proportional to the two-way attenuation loss so that the received signal 

20 would have a more or less desired (e.g., equal) signal over the received frequency spectrum. The 
controller 20 is further configured to sum the received signals and to apply a matched filter to the 
signals to recover structural interfaces. The same can also be performed in a single transmitted 
ultrasound burst where the frequency is scanned over a desired range while modulating the 
amplitude in a desired manner to compensate for attenuation differences associated with the 

25 different frequencies. Frequencies between or within pulses can vary in virtually any manner, 
e.g., linearly increasing or decreasing, non-linearly increasing or decreasing, monotonically 
increasing or decreasing, randomly, etc. The controller 20 is further configured to control the 
delay circuits 26 to vary the phase between or within pulses. The controller 20 can cause the 
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delay circuits 26 to cause different pulses to have different phases, and/or single pulses to have 
different phasing within the pulse. 

Power levels for diagnostic applications vary depending upon the application. Typically, 
however, for diagnostic applications the burst duration is 1-100 cycles with a frequency between 
5 about 0.1MHz and about 5MHz, and preferably between about 0.5MHz and about 3MHz, 
although these values are exemplary only, and not limiting as to the scope of the invention. 

The reflected signals received by the transducers 16 are converted to electric signals for 
processing by the imager 12. The imager 12 is configured to process the electric signals to form 
an image of the scanned region. The images can be provided to the controller 20, and may be 
10 obtained simultaneously with other uses of the array 14 (e.g., by intermittently applying energy 
for obtaining images, or by using some of the elements 16 for imaging and other elements 16 for 
other purposes). 

The controller 20 is configured to analyze data from the array 14 and images from the 
imager 12. The controller 20 can process the reflected signals, analyzing the returned signals 

15 (or lack thereof) received by the elements 16 to identify indicia characteristic of a particular 
condition, e.g., fluid in a patient's sinus cavity. 

FIGS. 3-5, with reference to FIG. 2, show various configurations of the system 10, with 
only a portion of the system 10 shown, of particular use in diagnostic applications. For example, 
referring to FIGS. 3-4, a single transducer or a small array of transducers is used as an ultrasound 

20 transmitter 32 in lieu of the array 14. The configuration shown in FIG. 3 is preferably used for 
diagnostic applications where an image is not desired, while the configuration of FIG. 4 lends 
itself for use with or without producing images for analysis. For example, as shown, ultrasound 
shear waves may be directed through a material 34, such as a bone, in which there is a cavity 36 
(such as a sinus cavity, or an abscess, etc.). Reflections may be analyzed by the controller 20 to 

25 determine various properties, such as whether the cavity 34 is fluid filled. FIG. 5 shows a 
configuration similar to that shown in FIG. 4, with electronic or mechanical scanning of 
transmitted shear waves 38 indicated by an arrow 40. While the shear waves 38 are shown 
unfocused, the phasing of elements in the transmitter 32 could be adapted to focus the waves 38. 
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In operation, referring to FIG. 10, with further reference to FIGS. 2-9, a process 1 10 for 
performing diagnostics on, and/or providing therapy to, an object using the system 10 includes 
the stages shown. For exemplary purposes, the object is assumed to be the skull 28 as shown in 
FIG. 2, although this is exemplary only and not limiting to the scope of the invention. The 
5 process 1 10, is exemplary only and not limiting, and may be altered, e.g., by having stages added, 
removed, or rearranged. 

At stage 112, the array 14 is positioned relative to the object 28. The array 14 (or single 
element transceiver) may be positioned manually and/or by actuation of the positioner 23 by the 
controller 20. The array 14 is positioned such that it can transmit ultrasound energy to the 

10 surface of the object 28 to produce primarily shear waves in the object 28. The desired 

positioning may be determined in accordance with the angle of the surface of the object 28 as 
determined visually by an operator of the system 10. The angle of the surface of the object 28 
may also be determined by actuating the array 14 and imaging the surface of the object 28. The 
image of the surface may be used to reposition the array 14 and/or to electronically steer the 

15 beam(s) from the array 14 as appropriate. 

At stage 1 14, the transducer elements 16 are excited to produce a desired ultrasound beam 
directed at the object 28. The controller 20 selects the desired elements 16 for transmitting 
energy to provide the desired direction, coverage, and energy of the ultrasound energy to provide 
sufficient energy to a desired region or point to achieve desired results (e.g., for diagnosis, 

20 imaging, and/or therapy). The controller 20 regulates the frequency, phase, amplitude, and delay 
of the energy provided by the selected elements 16. The delay from element to element may vary 
to provide an electrically-controlled beam direction such that desired amounts of energy from the 
array 14 as a whole, and/or individual elements 16, are incident upon the object 28 at an incident 
angle between the longitudinal and shear critical angles, and are directed to reach a desired 

25 region in the object 28. The ultrasound energy is preferably provided in a coded manner as 
discussed above, with the energy being pulsed over a range of frequencies, with power 
magnitudes that vary in accordance with (depending on) the frequency of the signal transmitted. 
The beam may be broad or focused. Elements may also be selected and excited to provide 
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ultrasound waves that will produce primarily longitudinal, not shear, waves in the object. The 
beam transmitted by the array 14 is incident upon the skull as longitudinal waves, passes through 
the skull as shear waves, and is converted back to longitudinal waves for further transmission 
(e.g., in the brain). 

5 Stage 1 12 and/or 1 14 may be repeated depending upon the application. For example, for 

imaging, the array 14 may be repositioned mechanically and/or the direction of the ultrasound 
beam electronically steered to cover the entire area to be imaged. 

At stage 116, for diagnostic or imaging applications, the imager 12 and controller 20 
processes received reflections to determine appropriate characteristics. The imager 12 may 

10 process the received signals by collecting them in conjunction with the corresponding direction 
of the incident beam to produce an image of the object 28 (e.g., at a plane within the object 28). 
The controller 20 may manipulate the returned indicia to determine one or more properties not 
necessarily related to imaging, such as a determination that a cavity within the object 28 contains 
fluid or not. For diagnostic applications, the image(s) provided by the imager 12 may not be used 

15 by the controller 20, and indeed the imager 12 may not provide any images. 



Therapeutic Applications 

The system 1 0 may be used for therapeutic applications, in addition to or instead of 
diagnostic applications separately or simultaneously. The components of the system 10 generally 

20 operate as described above. For therapeutic applications, however, another imager 13 may be 
provided and the system 10 configured to provide ultrasound in a manner tailored for therapeutic 
use. The elements 16 are preferably configured to pulse, and direct, ultrasound with high power 
levels for therapeutic applications (the levels used depending upon the application). 

The controller 20 is configured to control the phase shifters 26 and the amplifiers 24 for 

25 therapeutic applications. The controller 20 can cause the transmitted ultrasound to be focused on 
a focal region if appropriate for the therapy (e.g., thermal ablation) or to be applied to a broader 
region in an unfocused manner for appropriate therapy (e.g., low-power applications such as 
opening the BBB). 
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Power levels for therapeutic applications vary depending upon the application. Typically, 
for therapeutic applications, the burst duration is 1 cycle to continuous wave with a frequency 
between about 0.1MHz and about 3MHz, and preferably between about 0.2MHz and about 
2MHz, although these values are exemplary only, and not limiting as to the scope of the 
5 invention. The sonications typically will include one burst or multiple bursts resulting in 
sonication times between about 1 microsecond and about 1 h or longer with varying times 
between pulses depending on the application. 

The imager 13 can obtain images of the object 28 independently of the other apparatus of 
FIG. 2. For example, the imager 13 may be a computer tomography (CT) or magnetic resonance 

10 (MR) imager. The imager 13, like the imager 12, may provide images to the controller 20, e.g., 
for use in determining if adjustments are needed before or during therapeutic application of 
ultrasound using the array 14. The controller 20 is configured to determine from the images 
whether a desired region for receiving therapy is indeed receiving therapy, and receiving 
sufficient energy for the desired therapy. If a region other than the desired, target region is 

15 receiving the ultrasound energy (e.g., being ablated), or the desired region is receiving too much 
or too little energy, the controller 20 can alter the excitations of the elements 16 to correct the 
situation. For example, the energy can be increased or decreased, the focus of the energy made 
smaller or larger, the focal region can be changed or moved, etc. 

FIGS. 6-9 show various configurations adapted for use in therapeutic applications such as 

20 thermal ablation, mechanical destruction of tissue or tumors, BBB opening, gene therapy, 
targeted drug therapy, acceleration of chemical reactions, aiding for diffusion of chemicals, 
hyperthermia, applications increasing temperature of a region, any other ultrasound interaction 
that has therapeutic value, etc. FIG. 6 shows a single, curved transmitter 42, that may be a single 
element, or an array of elements. The transmitter 42 is configured to produce a focused 

25 ultrasound beam 44 that is transmitted through an object 46 (e.g., a skull) as shear waves, and 
that is focused at a focal point or region 48. Note that the ultrasound beam 44, as beams in other 
configurations discussed, is made of longitudinal waves between the transmitter 42 and the 
object 46 and shear waves in the object 46. If the makeup of the object 46 changes (e.g., from 
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bone to a cavity, be it filled or not, e.g., with fluid or other matter such as brain tissue), the 
transmitted beam may return to longitudinal waves. This is the case when the object 46 is a 
skull. FIG. 7 illustrates a configuration similar to that of FIG. 6, but with two transmitters 42, 43 
configured to focus their respective beams 44, 45 on the focal point/region 48. The transmitters 
5 42, 43 may be different portions of a single array, or physically different, e.g., different arrays, 
different single elements, or a single element and an array. FIG. 8 illustrates a configuration 
similar to FIG. 7, but with incident ultrasound waves 50 also used that produce longitudinal 
waves in the object 28. This configuration may be used, e.g., when power is desired in addition 
to the power that the transmitters 42, 43 can provide. The longitudinal waves may be provided 

10 by the array 14, which the transmitters 42, 43 may be part of as well. FIG. 9 illustrates the use of 
a linear transmitter 52 (e.g., a single element or an array of elements) displaced from the object 
46 and transmitting non-focused ultrasound energy of a sufficient level for a therapeutic 
application (e.g., opening the BBB, release or activation of chemicals or genes or other materials, 
etc.). As shown, the transmitter is mechanically positioned to provide a beam 54 that is normal 

15 to the transmitter 52 and between the longitudinal and shear critical angles with respect to the 
object 46. 

In operation, referring to FIG. 26, with further reference to FIGS. 2-9, a process 1 10 for 
providing therapy to an object using the system 10 includes the stages shown. For exemplary 
purposes, the object is assumed to be the skull 28 as shown in FIG. 2, although this is exemplary 
20 only and not limiting to the scope of the invention. The process 210, is exemplary only and not 
limiting, and may be altered, e.g., by having stages added, removed, or rearranged. Stages 212 
and 214 are similar to those discussed above. 

At stage 216, the desired region of the object 28 is imaged. The imager 12 and/or 13 
is(are) used to image the desired region, and preferably the surrounding volume such that the 
25 focus accuracy of the applied ultrasound can be evaluated. 

At stage 218, appropriate corrections are made to the incident ultrasound energy for 
therapeutic applications. If images produced by the imager 13 (or by the imager 12 if the array 
14 is suitable to both image and provide therapy) indicate that the applied ultrasound energy 
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should be adjusted, then the controller 20 acts accordingly. The images may be analyzed, e.g., to 
see if a desired region, and only the desired region, is being heated and determining what, if 
anything, is less than desirable about the applied ultrasound. For example, the ultrasound energy 
may be focused poorly, or in the wrong spot, or too narrowly, or of too low or too great of power, 
5 etc. The controller 20 can cause the beam to be modified by adjusting the phase, amplitude, 

delay, and or frequency of excitation in accordance with the analyzed images to adjust/correct the 
applied energy to yield a desired focus and amount of energy in the object 28 at the desired 
location. The controller 20 can also cause the positioner 23 to alter the orientation and 
relationship of the array 14 relative to the object 28 as desired. Adjustments are preferably made 
10 before application of full power for therapy (e.g., thermal ablation or other non-reversible 
procedures). For example, less than full power may be applied, a preliminary determination 
made from the images made as to whether and what adjustments should be made, adjustments 
made, and then full therapeutic power applied by the array 14. Adjustments may also be made 
during application of full therapeutic power based on progress of the therapy. 

15 

Experiments and Experimental Results 

A parallel plastic plate phantom was used to test the algorithm in order to provided an 
idealized isotropic case that could readily be verified by experiment. Relevant values for the 
plastic are summarized in Table I. The ultrasound source was a 1.5 MHz focused transducer 

20 with a diameter of 12 cm and a radius of curvature equal to 16 cm. The source function for the 
algorithm was acquired by projecting a laboratory-measured pressure field from a plane near the 
geometric focus backward to the source. This field was measured over a 30 X 30 area with a 
spatial resolution of 0.5 mm. The source function was propagated through the plates using an 
algorithm described numerically in Appendices A and B. The distance from the source to the 

25 measured plane was 121 mm and the distance from the plastic inner surface to the source was 61 
mm. 

Similar methodology was used to propagate through ex vivo human calvaria (brain cages). 
Following the laboratory measurements, water was used as the interfacing medium with the 
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inner and outer skull surface. The ultrasound source for the transskull measurement was a 0.74 
MHz focused transducer with a diameter of 8 cm and a radius of curvature equal to 15 cm. 

The numeric algorithm is implemented in Matlab®, using matrix-based operations for the 
layers. Operations were performed on a 1GHz AMD-Based PC. A typical projection of a 
complex 128 X 128 matrix through five layers took approximately 30 s to calculate. 



Table I 



Material 


Density 
(kg/m 3 ) 


Abs 
L 


Abs 
S 


c L 

(m/s) 


C s 
(m/s) 


Thickness 
(XI 0" 3 m) 


Plastic 1 


1187 


45 


50 


2185 


1330 


11.8 


Skull 0 


NA 












Skull 1 


2186 


85 


90 


2850 


1400 


5.29 


Water 


1000 


0 


0 


1486 


NA 


Variable 



LABORATORY MEASUREMENTS 
Skull Registration and Density 

Data for the simulation was obtained from a digitized human head profile obtained using 
CT images (Siemens, SOMATOM, AH82 Bone Kernel). Both the coordinates of the skull 
surfaces as well as the internal density variation are obtained from these images. Scans were 
taken at 1 mm intervals using a 200 mm x 200 mm field of view. A polycarbonate stereotaxic 
frame was attached around each sample to allow the skulls to be attached to the array and provide 
a reference for the mechanical positioning system and the CT images. The calculation was 
performed only in bone lying within the beamwidth of the section being considered. Information 
about the shape and structure of an individual calvarium was obtained from combining the 
images, that returned intensities proportional to material density. 

Coordinates of points along the inner and outer surfaces of the skull were identified on 
an image using a threshold filter that searched for the innermost and outermost densities >1.4 
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gm/cm along each line of an image. Points of successive images were combined to give three- 
dimensional representation of the inner and outer skull surface. Pixel intensities of each image 
were also combined into a three-dimensional array for later processing. 

The phasing algorithm relied on precise knowledge of the orientation of the skull relative 
5 to individual array elements. To achieve this task in practice, the phasing algorithm translated 
and rotated the skull data from the CT coordinate frame to the transducer coordinate frame as 
well as translated and rotated the skull from the mechanical positioning system's coordinate 
frame to the transducer coordinate frame. The program operated using three markers located on 
the polycarbonate frame affixed to the skull. These locations could be identified mechanically 
10 with the positioning system to a precision of approximately 0.1 mm. The algorithm generated a 
rotation matrix that mapped between the coordinate systems. 

Ultrasound Measurements 

Propagation experiments were set up in a water tank to verify the numeric algorithm. 

15 Measurements were performed in degassed and deionized water in a tank padded with rubber to 
inhibit reflections. Ultrasound signals were generated by a transducer specific to the particular 
measurement and received with a Polyvinylidene difluoride (PVDF) needle hydrophone 
(Precision Acoustics, Dorchester, UK). To help assure a strong reception while maintaining 
accuracy, a 0.2 mm-diameter hydrophone was used for 1.5 MHz measurements and a 0.5 mm- 

20 diameter hydrophone was used for measurements at 0.74. MHz. The smaller hydrophone at 
higher frequency was used to help minimize hydrophone directivity and to help prevent phase 
averaging. A skull (or plastic plate) was placed between the hydrophone and transducer at an 
angle controlled by a rotational stepping motor (made by Velmex of Bloomfield, New York). A 
3D linear positioning system (Velmex, Bloomfied, Model V P9000) allowed the hydrophone to 

25 be scanned over a measurement area centered about the transducer's axis of symmetry. 

Transducer signals were generated by an arbitrary waveform generator (made by Wavetek, of 
Norwich, UK, Model 305) fed to a power amplifier (made by ENI, of Rochester, New York, 
Model 2100L). The hydrophone's voltage response was sent though a Precision Acoustics pre- 
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amp and an amplifier (made by Preamble Instruments, of Beaverton, Oregon, Model 1820) 
before it was recorded by a digital oscilloscope (made by Textronix, of Wilsonville, Oregon, 
Model 380). The voltage waveform was downloaded to a PC by GPEB (General Purpose 
Interface Bus) control and the amplitude and phase at the measurement location were calculated 
5 from the FFT of the signal, taking values at the transducer driving frequency. 

The 1 1 .8 mm plastic plate was placed in the test tank and the acoustic transmission was 
measured on the transducer's axis of symmetry, 121 mm from its face. The acoustic pressure was 
measured between -70° and 70° at an increment of 1°. Agreement between the measured and 
simulated waveforms was evaluated by comparing the amplitudes and phases at each angular 

10 orientation. The source function for the simulation was a pressure field measurement taken with 
the transducer in water without a plate present. 

For the transskull measurements, two different mounting procedures were performed. 
The first, shown in FIG. 1 1 was designed to allow measurement through an approximately 
constant location on the skull at different incident angles. The second configuration, shown in 

15 FIG. 12, allowed measurement at high incident angles with good registration between the 
transducer and the skull, but did not facilitate movement of the skull. Initial measurements 
examined the angular dependence of amplitude on the skull, in order to determine whether 
amplitude peaks were present beyond the longitudinal Snell's critical angle. A section of skull 
bone skull was rotated between 0° and 55°, with the axis of rotation normal to the transducer axis 

20 and coincident with a line through the bone. The maximum angle of 55° was the highest value 
obtainable with the setup. 

RESULTS 

Plastic Phantom 

25 The correlation between the measured and predicted fields is shown in FIG. 13. For 

reference, the plot also includes the amplitude of the numeric calculation obtained when shear 
waves were neglected. The amplitudes and phases of all three curves match closely for incident 
angles below 31°, which is the longitudinal critical angle for the signal's spectral peak. However, 
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beyond this angle, the Longitudinal-only simulation is unable to predict the second local maxim 
in the amplitude, which results from a purely elastic wave within the sample. The major source 
of discrepancy between the measured data and simulation may be from result from 
underestimation of the shear wave absorption coefficient in the simulations. This effect becomes 
5 more pronounced at high angles, where the path length is longest. Very good correlation is 
found, however, between the simulated and measured ultrasound phase for all angles, with the 
exception of the transitional region between 20° and 30° where the wave amplitude is near its 
minimum. Neglecting this region, 76% of the remaining calculated points deviated from 
measurements by tc/6, or less, radians. 

10 

Transskull Measurement 

By rotating the skull, a local pressure transmission maximum was observed at 32°, which 
may result almost entirely from shear propagation through the bone, based on a longitudinal 
critical angle of about 30°. Further evidence of the wave's origin as a shear wave at higher 

15 incident angles was found in measured data shown in FIG. 14A-B, that show the angular 

dependence of the wave phase. As shown in FIG. 14B, below 26° the phase closely resembles 
the behavior of a purely longitudinal wave. Between 26° and 32° there is a deceleration in the 
slope, indicating a region of superimposed contribution from the shear and longitudinal wave. 
Beyond 32° the slope approximately linear and negative since the shear sound speed is less than 

20 the speed of sound in water. 

Using the same mounting configuration, field measurements were performed over a 30 
mm x 30 mm area with 1mm resolution, for different incident angles. FIG. 15 shows the field 
immediately after passing through the skull bone for incident angles of 0°, 15° and 35°, showing 
less distortion in the signal propagated through the skull as a shear wave than when propagated in 

25 a longitudinal mode at both 0° and 15°. 

A second set of measurements was performed with skulls placed in reference frames for 
spatial correlation with CT images. Although only an estimated value of the shear wave speed 
was used for the study, accurate spatial registration and data for the longitudinal modes allowed 
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the longitudinal critical angle to be identified. Skulls were aligned to assure that the transmitted 
signal was from shear propagation in the bone. FIG. 16 shows an example where the outer 
surface of the skull is oriented at 32° relative to the axis of symmetry, as determined from the 
simulation algorithm. Although there is discrepancy between the simulation and the 
5 measurement, the simulation is able to identify the presence of the focused shear wave. 

A demonstration was performed with a skull placed in a reference frame, to allow spatial 
correlation with CT images. Accurate spatial registration and data for the longitudinal modes 
allowed the longitudinal critical angle to be identified. Only an estimated value, however, of the 
shear wave speed was available for the study. Skulls were aligned to assure that the transmitted 
10 signal was from shear propagation in the bone. FIG. 7 shows a line measurement of the field, 
with the outer surface of the skull oriented at 32° relative to the axis of symmetry. Although 
there is discrepancy between the simulation and the measurement, the simulation was 
nonetheless able to identify the presence of the focused shear wave, with a beam shape 
resembling that of the measurement. 

15 

Diagnostic Applications - Transskull Imaging 

For transskull imaging, the reflected signal propagates a second time back through the 
skull, returning with a measurable signal, preferably with little distortion. Signals were analyzed 
from a 1MHz Panametrics imaging probe and pulse-receive amplifier. A nylon nut (c = 2.6 X 

20 10 6 mm/s Z = 2.9M Rayl) served as the imaging object. The nut was suspended above the floor 
of a water tank by a thin steel needle. A-line images were acquired with the propagation 
direction first perpendicular to the skull and then at about 33°. 

To help overcome losses, especially strong signal attenuation above 1 MHz, a new 
excitation scheme was used that sent a series of extended bursts through the skull. Bursts were 

25 repeated over a range of frequencies and the signals were numerically combined. In the 

exemplary experiment, 10 cycles were sent through the skull for frequencies ranging from 0.5 
MHz to 2 MHz at 0.01 MHz intervals. This series of measurements is shown in FIG. 17A. The 
signals were then normalized by dividing by peak value to give even weight, to accentuate the 
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weaker frequency signals as shown in FIG. 17B. Following acquisition, all of the time signals 
were summed and a matched filter was used to recover structural interfaces. The resulting signal 
is shown in FIG. 18 A. This figure demonstrates significant improvement in localization and 
signal amplitude as compared with a single burst. Having a normalized spectrum across the 
5 frequency domain should allow better localization of interfaces after matched filtering is applied. 
Furthermore, forced driven oscillations of the transducer likely will increase the transducer's 
bandwidth at the high and low ends of its response. Also, the combined numeric signal will have 
a net power far exceeding powers that could safely be achieved with a single burst. Thus, 
significant improvement in the signal-to-noise ratio is likely. 

10 The A-line images were followed by C-Scan imaging experiments using a focused 

broadband transducer array with a center frequency of 0.91 MHz. To construct the image, the 1 
MHz transducer was affixed to a stepping motor positioner. The diameter of the transducer was 
12.7 mm, representing the resolution limit of the present experiment. A scan of the nylon nut 
was obtained over a 20 mm x 20 mm area first in water, then through a skull fragment at 

15 approximately normal incidence, and finally at an angle of approximately 33 degrees. The 
combination of axial (temporal) and radial information by beam steering allowed three- 
dimensional images to be constructed and evaluated. Image slices of the object are shown in 
FIG. 19 along lines of constant time. Spatial resolution of the image was limited by the diameter 
of the present imaging transducer, which detected backscattered data directly above the element. 

20 A smaller beamwidth is expected to provide considerable improvement in image 

sharpness. Distortion between the water scan and the transcranial scan is at least in part caused 
by the method of image acquisition procedure. The act of moving the transducer over the skull 
surface introduced additional distortion due to the spatial variation of the skull thickness and 
acoustic properties. Aberration correction may be unnecessary in cases where a narrow 

25 ultrasound beam is directed through the skull. If a larger area array is used, however, 
noninvasive aberration correction algorithms could be employed. 

Noninvasive phase aberration correction could also be used to restore the focus through 
the skull bone. In this case, phase distortion could be predicted using carefully registered 
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information from CT scans of the ex vivo skulls. However, with the imaging hardware it is also 
possible that the overall phase distortion can be predicted from signature reflections off of the 
base of the skull. The reflection point would serve as a virtual transmitter in the brain and could 
be used for focusing at any point in the brain. Success with this method could improve focusing 
5 and potentially eliminate the need for CT scans of the head. In practice, a transducer array with a 
center frequency between about 0.5 MHz and about 2 MHz and with at least 100 elements may 
be used to help assure ample beam steering ability. 

Improvements in the signal could be quantified in terms of SNR as well as measurements 
of the reduction in object distortion and the accuracy of object locations. The apparatus setup 
10 (e.g., that shown in FIG. 2) preferably uses a fast switching multiplexer to allow waveforms to be 
recorded from multiple channels across the array. A comparison of the images at 0 degrees 
incidence and 33 degrees is presented in FIG. 20, showing a reduction in the image noise at 33 
degrees. 

15 Diagnostic Applications - Cavity Imaging 

Application of the transskull method may be used in applications not only for the brain, 
but also for the sinus cavities within the skull, that fill with fluid when infected. Referring to 
FIG. 21, to demonstrate applicability, a plastic phantom 130 was developed that provided a 
cavity 134 that could be filled with fluid (water) and drained. A 1-MHz center frequency 

20 Panametrics transducer 132 was pulsed in order to send a signal through the phantom 130 in one 
of four situations: (1) A high angle of ultrasound incidence (>45°) to induce shear wave 
propagation through the phantom 130 with an the cavity 134 being air filled; (2) A high angle of 
ultrasound incidence with a the cavity 134 being water filled; (3) A low angle of ultrasound 
incidence (<5°) giving nearly all longitudinal wave propagation through the phantom 130 with 

25 the cavity 134 being air filled; and (4) A low incidence angle with the cavity 134 being water 
filled. 

An example of advantages for cavity imaging using shear waves is presented in the 
results of the experiments shown in FIGS. 22A-B. FIG. 22A shows a dramatic difference in the 
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time history of the high angle signal when the cavity is air filled, and when it is fluid filled. In 
contrast, the longitudinal results shown in FIG. 22B are complicated in both the air and fluid 
filled cases by additional signal information coming from reflections within the bone itself, 
making accurate diagnosis of the imaging situation difficult. These results suggest that in 
5 practice, the shear mode could have clinical application as a method of determining the presence 
of fluid in sinuses. 

DISCUSSION 

Preliminary simulations and experimental measurements indicated that a coherent, 

10 focused ultrasound beam can be transmitted through the skull as a purely shear wave. At a 

driving frequency near 0.7 MHz, the peak amplitude through the skull due to shear propagation 
in the bone was found to be on the order of, and sometimes higher than, longitudinal 
propagation. Furthermore, the fact that the shear wave experienced a reduced overall phase shift 
confirms the viability of simplifying and extending noninvasive transskull focusing methods to a 

15 larger region in the brain. Although the shear amplitude was lower than longitudinal modes, 
there was evidence that the overall beam may suffer less distortion when propagated through a 
localized region of the bone. Furthermore, the shear mode experienced less phase distortion, 
presumably due to the similarity between the shear wave speed and the speed of sound in water. 
Comparable similarities would be found in soft tissues. More precise correlation between the 

20 simulation and measured data may be obtained using a more accurate measurement of the shear 
wave speed in the skull bone. 

There are several possible direct implications of intentionally using shear waves in 
transskull applications (although the invention is not limited to this application. First, shear 
wave propagation may be added to pre-existing noninvasive transskull phasing algorithms for 

25 improved focusing at high incident angles. This is particularly important when focusing close to 
the skull surface, where high angular incidence is used. Second, shear wave propagation may 
have application in transskull imaging, where a narrow ultrasound beam would be directed 
through the skull at intentionally high incident angles. Third, shear wave propagation could 
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potentially be used to detect flow in the brain by means of frequency Doppler shifts. All of these 
techniques may be performed in the sub-megahertz frequency range discussed. 

Other embodiments are within the scope and spirit of the appended claims. For example, 
due to the nature of software, functions described above can be implemented using software, 
hardware, firmware, hardwiring, or combinations of any of these. Features implementing 
functions may also be physically located at various positions, including being distributed such 
that portions of functions are implemented at different physical locations. Further, if different 
pulses are used having different frequencies and amplitudes, each of the pulses need not have 
only one frequency and one amplitude; one or more of the pulses may have more than one 
frequency and/or amplitude. Also, incident ultrasound energy need not be directed at a subject 
above the longitudinal critical angle in order to induce shear waves. The incident angle could be 
below Snell's angle, e.g., normal to the surface of the subject and/or the subject's bone. It is, 
however, preferred to direct the energy to be incident at an angle between the longitudinal and 
shear wave critical angles as this may help produce a desired amount of shear wave energy. 

APPENDIX A 

The pressure amplitudes were calculated by solving for the velocity potentials described in 
Eqns. (1), (2), (3) and ( 8). After linear algebraic simultaneous solution of the equations, it may 
be shown that the shear and longitudinal wave amplitudes in the bone are given by. 



A' s '(0') = 



-4A' L c' 2 k' 2 kj;p' cose 7 cos6 f sine/ 7 /_k[ 
D„ (0' ) p" cos0' + C// V p"C n (0' ) 



(9) 



and 



A , l{0 , ) = 
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where, 
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B„ (© ' ) = k'J cos 26/ 7 + 2k'[ sin9/' sinG / , (11) 
C„ (Q') = k s cosG/ 7 cos 2Q"+k L sin 26 77 sine /' , (12) 

D„ (0 ' ) = fc/ 7 ^ 7 (cf cos 26/' - 2c 'J 2 sin6 / 7 sin (e/ 7 - 26/ 7 )) , (13) 

and the transmitted angles are understood to be functions of 0, related by SnelPs Law, 

sin<9 7 sin#, 77 sine 77 
c c L c s 

The amplitudes of the waves transmitted into the brain from the incident longitudinal and 
shear waves are equal to 

Alu , = - {^s (6 >^P 5 ,cos9/% cos 29/' + sinefjsinief )/^ w 

Z)(0')p" cos0 /7/ +c ///2 ^ /// p /// C // (^ / ) 



and 



f _ 2 AH (9 ' 2 p " cos9, (cf 2 - cf 4- of cos 29/' )*„ (9 ' ) / 



10 A u L [{G l ) = ±— 1 (16) 

p u cos*'"-/),, (ff') + c ff/2 ty zosO'lB„{e') 



APPENDIX B 

Regardless of whether or not the layers are parallel, the transfer function may still be 
written in a closed form. For the thickness across the z-axis, z n , the sound speed c n , and density 
15 of each layer (n=I, n, EQ), that the unit vectors normal to the layer surfaces h n is calculated. For a 

given initial wavevector k 0xy , the ray path from (0,0, z 0 ) between any two surface interfaces 
traverses a distance of 

K nxy\- r — : > 
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where, as depicted in FIG. 2, F nxy is the vector extending along the layer from z-axis to the 
intercept of the layer with the ray. The unit vector along the wavevector' s path is given by k . 
Again, the frequency dependence on the wavevector orientation is understood. It follows that the 
ray position vector must be equal to R^, = R^ k )ixy . Although the initial wavevector orientation, 

5 k 0xy , is known, the direction of the wave vector in the first and subsequent layers must be 
calculated using the relation 

7 (*.♦, x ) = t-(" w+i * Luy ) > (18) 

which is a consequence of Snell's law in three dimensional space that requires the 
incident wavevector, the transmitted wavevector and the normal vector all be in the same plane. 
10 The transmitted wavevector on the right hand side of Eq. (18) may be obtained by crossing both 
sides of the equation with n n+l . Using cross product relations, it may be shown that 

Li* = 7 ±L (*- V - Ms y^') +C0S Tu^ ( 19 > 

The incident unit wavevector of the n th layer is equal to the transmitted wave of the (n-l) st 
layer. With the exception of the 0 th layer shown in FIG. 2, r nxy is found in order to calculate Eq. 

15 (17). Given the thickness across the z-axis of each layer, the point of intersection of ray R n with 
the surface of the n+1 layer is 

r*+\=R*+r H -**- (20) 

Over a series of N layers, the phase of a ray reaching the N th plane § R (k Nxy , 9 <o) is the sum 

of the phase contributions over each path length given by Eq.(l 7). The spatial phase at the plane 
20 z, is related to the ray phase at N by 

k N r N siny Nxy , as illustrated in FIG. 2. A ray leaving the 
initial plane with a polar angle y 0 ^, (co ) will arrive at the plane z with a new orientation y Nxy (a>) 
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determined by the N-l unit vector k t 




(21) 



given $(k Qxy9 <n,z Q ) , the phase of p(k x9 k y9 G> 9 z Q ) at the initial plane. 



The pressure over the plane at z can be expressed in terms of the ray phase presented in Eq. 
(2 1 ) and the transmission coefficient given in 



where the terms in square brackets on the right hand side of Eq. (22) can be viewed as an 
operator that maps the field from ko-space to a new k N -space. In the present problem, the 
amplitudes, T, are given by to Eqns. (15) and (16) and the component angles,-/ N . In practice, 

this mapping requires interpolation to produce a linearly spaced matrix at z. 
What is claimed is: 



p{k Nxy 9 G) 9 z) = p(k 0xy , co , z 0 )e 



(22) 



42 



